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Do binaries in clusters form in the same way as in the field? 
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ABSTRACT 

We examine the dynamical destruction of binary systems in star clusters of different 
densities. We find that at high densities (f 0'* — f 0^ Mq pc"^) almost all binaries with 
separations > 10"^ AU are destroyed after a few crossing times. At low densities 
(O(1O^)M0 pc^'^) many binaries with separations > fO'^ AU are destroyed, and no 
binaries with separations > 10^ AU survive after a few crossing times. Therefore the 
binary separations in clusters can be used as a tracer of the dynamical age and past 
density of a cluster. 

We argue that the central region of the Orion Nebula Cluster was ~ 100 times 
denser in the past with a half-mass radius of only 0.1 - 0.2 pc as (a) it is expanding, 
(b) it has very few binaries with separations > 10"^ AU, and (c) it is well-mixed and 
therefore dynamically old. 

We also examine the origin of the field binary population. Binaries with sepa- 
rations < 10^ AU are not significantly modified in any cluster, therefore at these 
separations the field reflects the sum of all star formation. Binaries with separations 
in the range 10^ - 10^ AU are progressively more and more heavily affected by dynam- 
ical disruption in increasingly dense clusters. If most star formation is clustered, these 
binaries must be over-produced relative to the field. Finally, no binary with a separa- 
tion > 10"* AU can survive in any cluster and so must be produced by isolated star 
formation, but only if all isolated star formation produces extremely wide binaries. 



Key words: binaries: general 
methods: A^-body simulations 



stars: formation - open clusters and associations 



1 INTRODUCTION 

■ Most stars i n the field are observed to be in multiple 
' syste mfl (e.g. iDuauennov fc Mavoilll99ll : iFischer fc Marcvl 
1 19921 ). The binary fraction of young stars ap- 
pears to be e ven higher than tha t of the field (e.g. 
Mathieul Il994l: iPatience et al.l |2002| : iGoodwin fc Kroupal 



20051: iGoodwin et al . 
Kroupa fc Burkertl (| 2001 



20071 . and references therein). 



show that in dense clusters it is 
very difficult to make significant numbers of binaries from 
initially single stars and that it is impossible to widen an 
initially narrow separation distribution to the observed 
wide distribution. Therefore, the field binary fractions and 
properties must in some way mirror the outcome of star 
formation. Thus, one of the most significant constraints 
on models of star formation is that these models correctly 
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^ iLadal ||2006| ) points out that the binary fraction of M-dwarfs is 
low and so most stellar systems are single. However, the fraction 
of stars (rather than systems) that are in binaries could be more 
than half. 



predict t he numbers and prop erties of binary and multiple 
systems (|Goodwin et al.ll2007h . 

However, at least 50 per cent (probably 70 - 90 per cent) 
of stars are thought to form in clusters (e.g. iLada fc Ladal 
[2003) . It is known that internal dynamical processes in clus- 
ters can rapidly and significantly a ffect t he proper ties of 
binari es in clusters (e.g. [Heggie 19751: lHillsiri975 a b: Kroupal 
Il995al lbt iKroupa. Petr fc McCaughreanI Il999l )". The degree 
to which binary properties will be altered depends also on 
the density and lifetime of the cluster in which the binaries 
are born. 

The field binary population is thus a mixture of binaries 
from different density environments (which may or may not 
be the same initial population) which have been dynam- 
ically processed in different ways. Therefore, the outcome 
of (hydrodynamical) star formation simulations should not, 
and cannot, be directly compared with the field population. 
This paper is the first in a series in which we attempt to 
reconstruct the initial binary population resulting from the 
early stages (Class O/I) of star formation. This paper will in- 
vestigate the dynamical processing of binaries in clusters of 
different densities and examine the extent to which binaries 
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of different separations are processed. This will allow us to 
discuss the main aspects of attempting to model the origin 
of the field, and reconstruct the initial binary population. 

In this paper, we investigate how a field-like initial bi- 
nary population is modified in clusters. In Section [2] we de- 
scribe our initial conditions; we present our results and dis- 
cussion in Sections [3l|4] and [H and we conclude in Section |6l 



2 METHOD 

2.1 Initial conditions 

The simulated star clusters have masses ~ 10^ - 10^ M© 
and a range of half-mass radii of 0.1, 0.2, 0.4 and 0.8 pc. 
This means we can simulate clusters with a range of densi- 
ties from ~ 50 Mq pc"'' to ^ 10^ M0 pc~"^ covering almost 
the complete range of probable initial cluster densities. We 
summarise the properties of the clusters that we simulate in 
Tabled] 

For each set of initial conditions we run an ensemble of 
at least 10 simulations which are identical apart from the 
random number seed used to initialise the positions, masses 
and binary properties. 

Our clusters are set -up as initially virialised Plummer 
spheres (iPlumnieil IQlll) accord ing to the prescription of 
lAarseth. Henon fc WielenI l| 19741 ) . The Plummer sphere pro- 
vides the positions and velocities of the centres of mass of 
systems - which may be single or binary systems. 

To create a stellar syst em, the mass o f the primary star 
is chosen randomly from a iKroupal l|2002l ) IMF of the form 



N{M) oc 



M 
M 



-^■^ mo < M/Mq < mi , 
"^■^ mi < M/Mq < m2 , 



(1) 



where mo = 0.1 M0, mi — 0.5 M©, and m2 = 50 Mq. For 
simplicity we do not include brown dwarfs (BDs) in our sim- 
ulations. The effect of dynamical processing on BDs in clus- 
ters will be studied in a future paper (for the results of 
recent observations of BDs in Orion see iMaxted et al. l|2008l : 
and for existing t heoretical work see iKroupa et al.ll2003l and 
iThies fc Kroupa"2008l). 

We then assign a secondary to the system depending on 
the binary fraction associated with the primary mass. 

For a field-like binary fraction we divide primaries 
into four groups. Primary masses in the range 0.08 
M/M p < 0.47 are M-dwar fs, with a binary fraction of 
0.42 jPischer fc Marcvl Il992l fl K-dwarfs have masses in 
the r ange 0.47 ^ M/M q < 0.84 and binary fraction of 
0.45 (iMavor et al.l Il992l ) and G-dwarfs have masses from 
0.84 M/Mp) <..^-^ ' ^'t^ binary fraction of 0.57 
l|Duauennov fc Mayor! Il99l'). All stars more massive than 
1.2 Mq are grouped together and assigned a binary fraction 
of unity, as massive stars have a much la r ger binary fraction 
than low-mass stars (e.g. Abt et al. '1990': Mason et al."l99d; 
iKouwenhoven et al.l 12005. 2007 : .Pfalzner fc Olczak .2007i. 
and references therein). 



^ We note th at the M-dwarf binary fr action is highly uncertain. 
The results of lFischer fc Marcvl lll992l ') are probably only appro- 
priate for stars with mass > 0.3 Mq, and the binarity and sepa- 
ration dis tributions b elow thi s mass may be quit e different (see 
especially rLadall2006l and also lMaxted et al.|[20o3) . 



Table 1. A summary of the different cluster properties adopted 
for the simulations. The values in the columns are: the number 
of stars in each cluster (A^stars), the typical mass of this cluster 
(Mclustor)i the initial half-mass radius of the cluster (»'i/2)i the 
crossing time of the cluster (tcross ) and the initial binary fraction 
in the cluster (/bin)- 



^ stars 




ri/2 


Across 


/bin 


1500 


~ 10^ 


Mq 


0.1 pc 


~ 0.02 Myr 


100 per cent 


1500 


~ 10^ 


Mq 


0.1 pc 


~ 0.02 Myr 


field-like 


1500 


~ 10^ 


Mq 


0.2pc 


~ 0.05 Myr 


100 per cent 


1500 


~ 10^ 


Mq 


0.2 pc 


~ 0.05 Myr 


field-like 


1500 


~ 10^ 


Mq 


0.4 pc 


~ 0.1 Myr 


field-like 


1500 


~ 10^ 


Mq 


0.8 pc 


~ 0.4 Myr 


100 per cent 


1500 


~ 10^ 


Mq 


0.8 pc 


~ 0.4 Myr 


field-like 


100 


~ 10^ 


Mq 


0.4pc 


~ 0.4 Myr 


100 per cent 


100 


~ 10^ 


Mq 


0.4 pc 


~ 0.4 Myr 


field-like 



Clusters with an initial binary fraction of unity for all 
stars are also created in order to test the hypothesis that all 
stars form in binary systems and that single stars a re solely 
the re su lt of the dynamical proces sing of binaries (|Kroupal 
Il995al lbl: iGoodwin fc Kroupa|[2005l ). 



2.2 Binary properties 

Secondary masses are drawn from a fiat mass ratio distri- 
bution with the constraint that if the companion mass is 
< 0.1 Mq it is reselected. This maintains the underlying bi- 
nary fraction, but biases the mass ratios of low -mass systems 
towards unity (see IKouwenhoven et al.ll2009al ). We note that 
this means that we do not recover a Kroupa IMF despite our 
primaries being drawn from this distribution. 

Eccentricities of binary stars are drawn fro m a thermal 
eccentricity distribution (|HeggiQil975 : lKroupallT995a . 20081) 
of the form 

Me) = 2e. (2) 

The generating function for or bital periods are the 
l og-no rmal distributions obs e rved by iDuguennov fc Mavorl 
(|l99ll ) and lFischer fc Marcvl (| 19921 ) of the form 



/ (logP) = Cexp 



-(logP-logP) 



(3) 



where logP = 4.8, (Tiogp = 2.3 and P is in days. The pe- 
riods are then converted to semi-major axes. Binary sys- 
tems with small semi-major axes but large eccentricities 
ar e expected to undergo t idal c ircularisation, as observed 
in iDuquennov fc Mavoij 's (|l99ll ) sample. We include the 
effe cts of ti da l circu larisation in our simulations by utilis- 
ing [KrouEi's (Il995bl ) ' Eigenevolution' process. However, we 
note that the systems that are circularised by eigenevolution 
are so hard as to be generally unaffected by dynamical inter- 
actions and so the details of eigenevolution are unimportant 
(we also ran simulations without any eigenevolution and also 
with a simplified circularisation mechanism and found no 
significant differences) . 

By combining the primary and secondary masses of the 
binaries with their semi-major axes and eccentricities, the 
relative velocity and radial components of the stars in each 
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system are determined. These are then placed at the cen- 
tre of mass and centre of velocity for each system in the 
Plummer sphere. 

Sim ulations are run using the kira int egrator in Star- 
lab (e.g. iPortegies Zwart et al.l Il999l . l200ll . and references 
therein) and evolved for 10 Myr. Note that we do not include 
stellar evolution; we have checked that this makes very lit- 
tle difference to the outcome of the simulations as there are 
few high-mass stars that evolve on such a short timescale 
(and, as we shall see, most binary evolution occurs in the 
first Myr, before the massive stars have evolved). 



2.3 Finding bound binary systems 

We determine whether a star is in a binary system using 
the following methoc(f|. For each star, the distances to its 
ten nearest neighbours are determined. These distances are 
then used to find the average distance between stars locally. 
We then determine the identities of the nearest and sec- 
ond nearest neighbours of each star. If two stars are mutual 
nearest neighbours with a separation of less than a quarter 
of the local average distance between stars they are a po- 
tential binarjQ. Numerical experiments have shown that it 
is extremely rare to find a bound system with a separation 
greater than a quarter of the local average separation. They 
are a true binary if they also have a significantly negative 
relative energy, and they are an observational binary if they 
have positive energy (i.e. an observer may think the system 
is a binary when it is just two stars passing). Clearly in our 
definition an observational binary can be found in either 2D 
or 3D, however only 2D observational binaries are of interest 
as they represent potential mistakes in actual observations. 
However, for the rest of this paper we ignore observational 
binaries as they are found to be extremely rare. 



3 CLUSTERS OF DIFFERENT DENSITIES 

In this section we will examine the evolution of the binary 
populations in clusters of different densities. 

Firstly, we will consider dense clusters with half-mass 
radu of 0.1 - 0.2 pc and masses of ~ 750 M© (iVstars ~ 1500) 
with densities of lO" - 10^ Mq pc"^. Then we will consider 
low-density clusters with large half-mass radii, or low num- 
bers of stars, with densities of 0(10^) M0 pc~^. 

These two density regimes cover the whole range of ini- 
tial cluster densities in probably equal proportions. If the 
cluster mass function is N{M) oc M"^ (|Lada fc Ladall2003l ') 
then each equal logarithmic mass bin contributes the same 
mass of stars to the field. Taking the range of cluster masses 
to be 10^ - 10^ M0, then an equal mass of stars is produced 
by clusters larger than Orion (~ 10^ M©) to the mass pro- 
duced by clusters smaller than Orion. The higher densities 
are probably fairly representative of larger clusters, and the 



^ Our algorithm almost exactly repr oduces the results of the inde- 
pendent binary finding algorithm of iKouwenhoven et al.l ll2009bl . 
in prep.). 

* This method is also able to find higher-order systems by ex- 
amining the relationships between sets of three or more stars. 
For example a triple is a system in which three stars are mutual 
nearest and second nearest neighbours. 



lower densities of smaller clusters (we will return to this in 
Section |4] where we examine the initial conditions of Orion). 

We might well expect the evolution of the binary pop- 
ulations in different density clusters to be very different as 
the position of the hard-soft boundary, flh-s, will depend 
on the number of systems A^s va in, and the half-mass radius 
ri/2 



, of the cluster f following iBinnev fc Tremaine|[l987l ). 



ri/2 

iVsvs 



(4) 



Thus the hard-soft boundary will shift by a factor of nearly 
100 between the least and most dense clusters. 



3.1 Dense clusters 

3.1.1 The evolution of the binary fraction 



In Figs. 1(a) and 1(b) we show typical examples of the evo- 



lution of the binary fractions with time for clusters with 
initially field and 100 per cent binary fractions, respectively, 
with half-mass radii of 0.1 pc. The two lines represent the 
evolution of the binary fraction of all stars (dashed line), 
and of M-dwarfs (solid line). We define the binary fraction 
to be the fraction of multiple (almost always binary) systems 
compared to the total number of systems. 



/bin 



B + T 



S + B + T 



(5) 



where S, B, and T (etc.) are the numbers of single, binary 
and triple (etc.) systems, respectively. 

For both distributions there is a rapid disruption of sys- 
tems in the first 0.1 Myr, in which the soft systems are 
destroyed. Once the least bound (i.e. wide) systems are de- 
stroyed, the clusters reach an equilibrium in which the bi- 
nary fractions remain roughly constant. Thi s is the same 
behav iour as found bv lKroupal (|l995al lb[). and lKroupa et all 
(119991 ). 

The destruction of binary systems occurs in roughly 
a crossing time, and hence the binary population is set by 
the denses t configuration of a cluster before any expansion 
can occur l|Kroupall2o6ol ). As long as a cluster spends a few 
initial crossing times in a dense configuration, that is enough 
to process the binary population. The current density of 
a cluster (as long as it is lower) will merely maintain the 
population created in the dense state. 

What is also clear in Fig. [1] is that there is a mass- 
dependence in the destruction of systems. Low-mass binary 
systems are more susceptible to destruction as they have 
a lower binding energy. The binary fraction of M-dwarfs is 
~ 5 per cent lower than the binary fraction of the other. 



more massive stars, despite the scenario plotted in Fig. 1(b) 
in which the various spectral types have the same initial 
binary fraction. 

The wide field-like log-normal distribution used to gen- 
erate the initial binary population is not recovered even 
before any dynamical evolution. This is first apparent in 
Fig. |l(b)[ in which we only obtain a binary fraction of 
~ 75 per cent rather than 100 per cent at time zero, be- 
fore any evolution takes place. 

It is important to note that the difference is not due 
to a flaw in our binary finding algorithm. Rather, the aver- 
age distance between stars in the cluster is around 2000 AU 
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Figure 1. The evolution of binary fractions with time for a 0.1 pc half-mass radius cluster with (a) an initial field binary fraction; and 
(b) an initial binary fraction of 100 per cent, in which all stars are formed in binary systems. The binary fraction for all stars is shown by 
the dashed (black) line and the binary fraction for M-dwarfs is shown by the solid (red) line. The initial binary fractions (at cluster birth) 
are denoted by the plus sign (all stars) and the asterisk (M-dwarfs). In both scenarios, a large proportion of binaries are immediately 
broken up before an equilibrium state is reached after a few Myrs. The final M-dwarf binary fractions after the cluster has evolved for 
lOMyrs are 25 per cent for an initial field binary fraction, and 50 per cent for an initial binary fraction of 100 per cent. 



within the half-mass radius. Therefore binaries with separa- 
tions of this order or larger - whilst generated by the initial 
conditions - are not physically associated or bound within 
a cluster, and therefore are not identified as binaries. 



3.1.2 An initially field-like binary fraction 

In Fig. [2] we show the binary fraction-separation distribution 
with initially field properties in a cluster with half-mass radii 
of 0.1 and 0.2 pc. The open histograms show the distribu- 
tion at time zero - before any dynamical evolution of the 
systems, the hash ed histograms those at 1 Myr (roughly 
the age of Oriom JeffriesI 2007a 3)- The s olid (red) lines 
are the G-dwarf (iDuguennov fc Mavoij[l99ll ) log-normal pe- 
riod distribution, and d ashed (blue) lines are the M-dwarf 
IIFischer fc Marcvlll992l ) log-normal distribution, which are 
the functions used to generate the initial binary populations. 
The open circles in Fig. |2(a)"| show the binary fractions gen- 
erated by the initial conditions (there is a deviation from the 
generating function at low separations due to the effect of 
eigenevolution, however, as we shall see later this is unim- 
portant as these binaries are hard). The open histograms 
show the distribution of binary fractions found by our bi- 
nary finding algorithm, which are clearly different. 

In Fig.[2]we also compare the separation distributions of 
binaries generated with a field binary fraction found by our 
binary finder at time zero (open histogram), and (hatched 
histogram) at 1 Myr. Quite clearly there has been signifi- 
cant dynamical destruction of binaries with separations of 
~ 100 - 1000 AU, whilst binaries with separations < 50 - 



ary in our clusters is a t a few hundred AU (t his is also seen 
by e.g. iKroupa 1995a b |: iKroupa et al.lll999l and described 
in detail bv lKroupall2008l )! 



100 AU are almost uncha nged. This is the Heggie-Hills law 
(|Heggielll975l : lHillslll975al lbi) in action: the hard-soft bound- 



Even the binary fraction of very hard systems is re- 
duced by the destruction of wide binaries. There is very 
little evolution in the number of systems with separations 
below 1 AU; however the binary fraction of those systems 
decreases due to the increase in the total number of sys- 
tems due to the destruction of wider binaries. For example, 
in a cluster with 100 binary systems, 20 of these may be 
very hard. The initial binary fraction of these hard systems 
would be 20/100 = 20%. However, after the destruction of 20 
wider systems, each wider system becomes 40 single stars, 
and so the binary fraction of very hard systems would be 
20/120 = 17%, despite none of them having been destroyed. 

There is very little dynamical processing of initially 
hard systems. In each cluster a few (0 - 5) systems with 
initial separations <^ 50 AU are significantly altered or de- 
stroyed, but most systems retain virtually unchanged sepa- 
rations from formation. 

We conclude, in common with other authors, that bi- 
naries in dense clusters cannot form with the field separa- 
tion distribution as many field binaries are too wide to have 
formed in a dense environment. A cluster in which the aver- 
age separation between systems is a few thousand AU cannot 
possi bly form systems with separations grea ter than this. In- 
deed, IScaUyj^^larke^^&l^cCaughre^ l|l999t ) find only 3 pos- 
sible binaries in Orion with separations of 1000 - 5000 AU 
(they also note that the origin of wide binaries cannot be in 
Orion-like clusters). 
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Figure 2. The evolution of the separation distribution for a cluster with an initially field binary fraction and half-mass radius of (a) 
0.1 pc; and (b) 0.2 pc. The separation distribution observe d for field G-dwarf s bylP ugucnnoy & Mavor (1991) is shown by the solid (red) 
log- normal; the distribution observed for field iVI-dwarfs bv lFischer fc Marcvl lll992l ) is shown by the dashed (blue) log- normal. The open 
circles in (a) show the initial distribution generated by our initial conditions generator, whereas in both panels the open histograms show 
the initial binaries and the hatched histogram shows the binaries remaining after 1 Myr, as found by our algorithm. 




Figure 3. The evolution of the separation distributions for clusters containing ~ 1500 stars created with a 100 per cent binary fraction. 
The open histograms are the initial distribution and the hatched histograms are the distributions after 1 Myr. We show the separation 
distributions fo r such clusters with initial h alf-mass radii of (a) 0.1 pc; and (b) 0.2 p c. The log-normal fits ob tained from observations of 
field G-dwarfs llDuquennov &: Mavor|[l99ll . the solid (red) line) and field M-dwarfs llFischer &: Marcvlll992l , the dashed (blue) line) are 
also plotted. 



3.1.3 An initially 100 per cent binary fraction 



Fig. 1(b) appears to show that if the initial binary fraction 
is unity in dense clusters, then the effect of dynamical evolu- 
tion is to lower the binary fraction to close to the field values 
(actually slightly too high for M-dwarfs). This might suggest 



that in dense clusters stars form with a field-like separation 
distribution, but with a higher binary fraction (e.g. unity). 

However, as we show in Fig. |3](c.f. Fig. [2]), exactly the 
same effects occur with a binary fraction of unity as with 
a field binary fraction. Firstly, many of our generated bina- 
ries are unphysically wide, given the cluster's size, and are 
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(a) 0.8 pc half-mass radius, 1500 stars, field binary fraction. (b) 0.8 pc half-mass radius, 1500 stars, 100 per cent binary frac- 
tion. 





Separation (au) Separation (au) 

(c) 0.4 pc half-mass radius, 100 stars, field binary fraction. (d) 0.4 pc half-mass radius, 100 stars, 100 per cent binary fraction. 



Figure 4. The evolution of the separation distributions for clusters with initial half-mass radii of 0.8 pc, 1500 stars, with field-like and 
100 per cent binary fractions ((a) and (b) respectively); and clusters of initial half-mass radii of 0.4pc containing ~ 100 stars created with 
field-like and 100 per cent binary fractions ((c) and (d) respectively). The open histograms are the initial distrib ution and the hatched 
histograms are the distributions after 1 Myr. The log-normal fits ob tained from observations of field G-dwarfs ijDuauennov fc Mavoil 
I1991I . the solid (red) line) and field M-dwarfs llFischer fc Marcvll 19921 . the dashed (blue) line) are also plotted. 



not identified as binaries even before the start of the simu- 
lations. Secondly, the hard-soft boundary is in exactly the 
same place and so many binaries with separations > 50 to 
a few hundred AU are dynamically disrupted. 

Our initial conditions also produce too many low- 
separation binaries which are unaffected by dynamical evo- 
lution. This suggests that the initial hard binary population 
should look similar to the field (at least after mixing with 
clusters of lower densities and field stars - see Section [Sjl . 

It should be noted that high-density clusters - even 
though they start in virial equilibrium - expand significantly 



due to the input of energy into the cluster from the disrup- 
tion of significant numbers of binary systems. 



3.2 Low-density clusters 

We now examine the dynamical evolution of binary systems 
in low-density environments. In Fig.|4]we show the evolution 
of an initially field-like initial binary population in a cluster 
with ~ 1500 members with a half-mass radius of 0.8 pc, 
and of a cluster with ~ 100 members and a half-mass radius 
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of 0.4 pc. These clusters both have densities of ~ 200 - 
300 Mq pc~^ (compare with the clusters illustrated in Fig. [2] 
which have densities of ~ 10* pc~^). We show the binary 
separation distributions initially (zero time, rather than the 
initial conditions), and at 1 Myr for the clusters. 

What is clear in all panels in Fig. 3] is that most very 
wide binaries (> 10'* AU) never survive (even if they form) 
even in low-density clusters. The hard-soft boundary has 
clearly shifted in the low-density clusters compared to the 
high-density clusters considered above, but not enough to 
allow these extremely wide binaries to survive. 

From Figs. [2] and [3] (and inspection of the data) we 
estimate the hard-soft boundary in a dense A'' ~ 1500, 
fi/2 = 0.1 pc cluster to be ~ 50 - 100 AU. For the low 
density cluster with N ~ 1500, ri/2 ~ 0.8 pc, Eqn.|4] would 
suggest that the hard-soft boundary should be at around 
500 AU. However, inspection of Fig. |4] seems to show that 
the limit of the destruction of intermediate binaries is around 
5000 AU, about ten times larger than expected. 

The reason for this discrepancy is that the low den- 
sity cluster is still dynamically young at the 1 Myr ages 
illustrated in the previous figures. The crossing times of the 
dense clusters are only 10s kyr meaning that the clusters 
are many crossing times old. However, the crossing times of 
the low density clusters are around 0.4 Myr, meaning that 
they are dynamically very young. As illustrated in Fig. [T] 
and described above, it takes a few crossing times to reach 
an equilibrium in which all of the soft systems have been 
destroyed and the low density clusters have not had long 
enough to do this by 1 Myr. This also means that the low 
density clusters are not dynamically mixed and the binary 
fraction has a radial dependence (we will describe this in 
detail in the next section). 



4 THE INITIAL CONDITIONS OF THE 
ORION NEBULA CLUSTER 

The Orion Nebula Cluster (ONC) is observed to have a half- 
mass radius of ~ 0.8 p c and around 2000 - 3000 members 
jMcC aughrcan & Sta uffeil 1 19941 : iHillenbrand fc HartmannI 
ll998i : lKohler et al.i200^ . This makes the ONC a low-density 
cluster following our definitions above. However, this is the 
current state of the ONC. As shown bv lBastian et al.| (|2008| ) 
many clusters seem to undergo an early dense phase and 
then rapidly expand (po ssibly due to gas expulsion, see 
ICoodwin fc Bastianll200a ). In this section we examine the 
binary properties of the ONC and its dynamical state to 
atte mpt to infer the initial conditions of the ONC. 

iReipurth etHI |20o3) carried out a survey of inter- 
mediate/wide (visual) binaries in the ONC. They looked 
at the separation distribution in the ra nge of ~ 68 A U to 
676 AU. As noted in Section [3T2] iScallv et al.1 l|l999l ) find 
only 3 possible binaries in O rion with separations greater 
than 1000 AU. Therefore, the lReipurth et al.l (|2007l ') sample 
covers most of the range of intermediate to wide binaries 
in the ONC (presumably there are also a small number of 
binaries in the range 700 - 1000 AU). 

The immediate and obvious conclusion to draw from 
this is that ONC-like clusters - even at low-density - can- 
not be the source of a significant number of binaries with 
separations > 1000 AU in the field. 



The iReipurth et al.l (|2007l ') data are shown in Fig. (5] 
they place their data into bins of width 0.2 log AU and 
the numbers of binaries in each bin are shown by (green) 
crosses. We also plot the log-normal distributions for G- 
dw arfs and M-dwa r fs. No te that most of the stars observed 
by IReipurth et"al] (|2007l ) are M-dwarfs, and so it is the 
(lower) M-dwarf distribution that should be compared to 
the data rather than the G-dwarf distribution. 

We over-plot the average of our ensemble of simulations 
for clusters with 100 per cent initial binary fractions and 



half-mass radii of 0.1 (Fig. 5(a) I and 0.2 pc (Fig. 5(b) I, and 
an initially field binary fraction and half-mass radii of 0.4 pc 



(Fig. 5(c) I and 0.8 pc (Fig. 5(d) I, shown by the histograms 



in each panel. Each histogram shows the average for that 
bin of 10 simulations, with the dot-dashed lines showing the 
standard deviations over the 10 realisations. 

From inspection of Fig.[S]it is clear that the observations 
of the ONC can be matched by both an initially high binary 
fraction in a high-density cluster, or an initially field-like 
distribution in a low-density cluster. 

The current half-mass radius of Orion is ~ 0.8 pc, a fac- 
tor of ~ 5 larger than that inferred from an initially high bi- 
nary population. Despite the fact that a field-like population 
can match the observations we strongly prefer the interpre- 
tation that Orion started with a far higher binary fraction 
in this separation range which has dynamically evolved into 
the current distribution. The reasons for this are threefold. 

• Firstly, the central region of Orion (0-5 pc) is thought 
to be expandi njfl, based on the ID velocity dispersion of 



pcrsion is ~ 4.3 kms * ( 


Kroupa. Aarseth & HurlevI 200ll: 


lOlczak, Pfalzner & Eckart 


2008|). which is too large for it to 



be in virial equilibrium (although we note that the effect of 
binaries in a cluster is to incre ase the observed velocity dis- 
persion; [Ko^^nhove^^^^eTjlriis 2008) - a 0.8 pc half-mass 
radius cluster with a mass of 1500 Mp has a virialised veloc - 
ity dispersion of ~ 2.8 kms^* fsee lKroupa et aDll999l . [200ll ). 
Clusters are expected to expand early in their lives due to 
the effects of gas expulsion (see iBastian fc GoodwinI l2006l : 
Goodwin & Bastian 2006). Indeed, a factor of 6 increase in 
radius is not unusu al for a low effective star formation ef- 
ficiency cluster (see lCoodwin fc Bastianll2006l . in particular 
then fig. 2). The current velocity dispersion of 4.3 kms~* im- 
plies an initial size of the cluster of < 0.2 p(|f|. Therefore, its 
initial configuration must have been denser than is currently 
observed - and, as we have seen, it is the initial, rather than 
current density that is all important for processing th e initial 
binary population. IScallv. Clarke fc McCaughreanI (|2005| ) 
suggest from their detailed study of the dynamics of the 



^ On larger scales (~20pc), the region ar ound the ONC may be 
undergoing cold collap se, as observed by iFeieelson et al.l fcOOSl ') 
and Flbbin et all l|2009l) . 

^ Whilst an initially super- virial cluster expands its velocity dis- 
persion does decrease. However, to have expanded to 0.8 pc with 
a velocity dispersion of 4.3 kms~^ in 1 Myr the young ONC 
must have been significantly super-virial, probabl y with an ef- 
fective star formation efficiency (eSFE) of ~ 0.3 llKroupa et al.l 
I2OOII : iGoodwin fc BastianI [20061) . Such a low eSFE would sug- 
gest that the velocity dispersion would not decrease significantly 
during the expansion as th e cluster would be globally unbound 
dCoodwin fc Bastianll2006l) . 
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Figure 5. Comparison with fig. 7 from lReipurth et al.l l[2003). Their data are shown by the (green) crosses, together with their uncertain- 
ties. Data from our simulations for a cluster age of 1 Myr and cluster membership of ~ 1500 stars (similar to the ONC) are shown by the 
open histograms (uncertainties shown by the da sh-dot lines and offset from the centre of each bin for clarity) . As in ou r Figs. [H [3] and [H 
the l og-normal distributions for field G-dwarfs llDuauennov fc Mayor] [l99ll . the solid (red) line), and field M-dwarfs jFischer &: Marcvl 
I1992I . the dashed (blue) line) are plotted. We show the results for three different initial cluster half-mass radii. Panels (a) and (b) show 
the results for an initial binary fraction of 100 per cent and half mass radii of 0.1 pc and 0.2 pc respectively. Panels (c) and (d) show the 
results for an initial field binary fraction and a half mass radii of 0.4 pc, and 0.8 pc respectively. 



Orion cluster that it may well have initially been 10 - 100 
times denser in the past; we suggest that it was around 100 
times denser, in broad agreement with them. 

Secondly, the binary population observed by 



iKohler et all (120061 ) and lReipurth et al.1 ll2007h in Orion ap- 
pears to be well-mixed. iKohler et alT(|2006D find very little 
difference in the binary fraction (more correctly, the com- 
panion star frequency, see their figs. 10 and 11) with radius. 



This is difiicult to explain if Orion started large, as the in- 
ner regions should have been more heavily processed than 
the outer regions as the current crossing time for the clus- 
ter is of order its age (see e.g. IKohle r 2004*) . However, if 
the cluster started out ~ 100 times denser than it is now, 
then it is significantly dynamically older than its current 
size would suggest (as w ould be the case for expanding clus- 
ters; [Bastia^^et^ ] |2008i ). Hence its binary population would 
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Figure 6. The binary fraction as a function of radius after 1 Myr of evolution for clusters containing 1500 stars but with different initial 
binary fractions and different initial densities: (a) with an initial half-mass radius of 0.1 pc and a binary fraction of 100 per cent, and (b) 
with an initial half mass radius of 0.8 pc and a field-like binary fraction. We plot the binary fraction of all stars in the cluste r against 
the di stance from the centre in half-mass radii. This enables a direct comparison with the ONC data in figs. 10 and 11 in iKohler et al.l 
1I2OO6I') . 



be well-mixed and evolved within a few hundred thousand 
years (as discussed above) and no very significant differences 
between the binary fractions in the inner and outer regions 
would be expected. In Fig. |6] we show the radial variation 
of the binary fraction at 1 Myr in a dense (ri/2 =0.1 pc) 
cluster with an initial binary fraction of 100 per cent and 
low-density (ri/2 = 0.8 pc) cluster with an initial field-like 
binary fraction. (Note that the dense cluster expands, so 
that after IMyr ri/2 = 0.2 pc, whereas the sparse cluster 
retains its initial half-mass radius.) Our simulations show a 
trend towards an increasing binary fraction with radius for 
the low-density (dynamically young) cluster, and no trend 
in t he high-density ( dynam ically old) cluster. 

iReipurth et"al] (|2007l ) do find that the ratio of wide-to- 
close binaries increases with radius within the inner pc, and 
is fiat beyond this (see their fig. 9; close binaries are 66 - 
225 AU, and wide binaries are 225 - 670 AU). However, 
easily within the errors the ratios are fiat beyond the inner 
0.5 pc, and it is only within the inner 0.5 pc that there are 
significantly fewer wide binaries than close ones. This is not 
at all unexpected. While the Orion cluster has been expand- 
ing from its proposed denser initial state, it will still undergo 
dynamical processing as wider binaries are always more sus- 
ceptible to destruction in the inner regions of a cluster than 
in the outer regions. If Orion was always at its current size, 
then the wide binary population in the outer regions is very 
close to the initial population, whilst in the inner regions it 



has been processed (the situation in Fig. 6(b) 1. However, if 



Orion was originally far denser, then the wide binary pop- 
ulation has been processed in every region of the cluster, 
only somewhat more so in the inner regions (the situation 
in Fig. 6(a) 1 and at no location in the cluster do we see the 
initial wide bi nary populati o n. Th is leads to our third point. 
• Finally, IScallv et all (|l999l ) find very little evidence 



of binaries with separations > 1000 AU in the ONC. If the 
ONC formed at low density, it must have formed with almost 
no binaries with separations > 1000 AU as a low density 
cluster could not have been able to destroy a population of 
binaries this wide in 1 Myr. But the lack of binaries with 
separations > 1000 AU can easily be explained if the ONC 
was far more dense in the past as it would have destroyed 
the majority of binaries with these separations. In addition, 
that we do not see binaries with separations > 1000 AU 
even in the very outer regions of the cluster suggests that 
the entire cluster has been dynamically processed. 

Whilst none of these arguments show conclusively that 
the ONC has expanded from an initially dense stage, we 
think the y are at least an ind ication. However, we would also 
note that I Allison et all (|2009, ApJL, submitted) can explain 
the formation of the Trapezium if the ONC was once ~ 100 
times denser. 



5 DISCUSSION 

Binary fractions and multiplicities are clea rly strong con- 
strain ts on theories of star formation (ICoodwin et al.l 
|2003). However, it is often very difficult to compare 
the outcomes of theory with observations. In particu- 
lar, the outcome of star formation theories cannot be 
directly compared to the field (as attemp t ed by e.g . 
iGoodwin. Whitworth fc Ward-Thompsonll20o4 lBatefl2009h 
as (a) there is significant dynamical evolution in dense clus- 
ters which will alter the star formation products, and (b) 
the field is the sum of many st ar forming regions and modes 
(e.g. iBrandner fc K6hleilll99d ). 

We have examined two density regimes - high-density 
clusters (10* — 10^ M© pc~^), and low-density clusters 
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(0(10^) M0 pc~'^). We have argued that these two groups 
bracket the vast majority of densities in cluster-forming star 
formation. We have seen that neither low- or high-density 
clusters allow binaries with separations > 10* AU to exist. 
And high-density clusters do not allow binaries with separa- 
tions > 10^ AU to exist. However, observations of the field 
show that such binaries apparently do exist in significant 
numbers. This raises a number of questions about the origin 
of the field population. 



5.1 Where do binaries of different separations 
come from? 

In the following discussion we will examine the origin of 
the field binary population from clustered and isolated star 
formation. For simplicity we will assume that 20 per cent 
of star formation occurs isolated (ISF) (actually the figure 
is probably so mewhat less than t his and could be as low 
as 10 per cent; iLada fc Lad"all2003l ). and 40 per cent each in 
low-density clusters (LDC) and high-density clusters (HDC) 
(i.e. low- and high-mass clusters with a /3 = —2 power-law 
cluster mass function). 

We will also ass ume that the G-dwarf field separation 
distribution (i.e. the iDuguennov fc Mavorlll99ll Gaussian) 
also holds for other stellar masses. In particular, we assume 
it holds for M-dwarfs which might not be a part icularly good 
assumption as (a) the lFischer fc Marcvl l|l992l ) data is only 
in rough agreement w ith the G-dwarf distribution, and (b) 
the lFischer fc Marcvl data p robably onl y holds for M-dwarfs 
with masses > O.SM© (e.g. iLadallBo^ . However, these as- 
sumptions at least allow us to discuss the pertinent points. 

We can divide binaries into four broad categories based 
on their separations. Binaries with separations < 50 AU are 
'always hard' - no density of cluster significantly changes the 
separation distribution in this range. Binaries with separa- 
tions of 50 - 1000 AU are 'sometimes hard' - high-density 
clusters can destroy some of this population, but low density 
clusters and isolated regions cannot. Binaries with separa- 
tions in the range 10^ - 10* AU are 'soft-intermediates' - 
high density clusters destroy such binaries, and low-density 
clusters destroy some. Binaries with separations > 10* AU 
are 'always soft' - any cluster will destroy such binaries (if 
they could even form in the first place). 

• Roughly 50 per cent of binaries are 'always hard'. 
Such systems cannot be destroyed by all but the most ex- 
treme cluster densities and so the field population must rep- 
resent the sum of binaries formed in HDC, LDC and ISF 
star formation. 

Thus, the binary fraction and separation distribution 
below 50 AU must be a fundamental outcome of the star 
formation process. That is, a combination of HDC, LDC 
and ISF must produce around 30 per cent of G-dwarf j3 and 
15 - 20 per cent of M-dwarfs with a companion < 50 AU 
and the combined separation distribution in this range must 
match the field. 

• Around 10 - 15 per cent of binaries in the field are 
'sometimes hard'. Therefore star formation would be ex- 
pected to slightly over-produce such systems as many in 



''' i.e. 50 per cent of the 60 per cent of G-dwarfs in binaries. 



HDCs will be destroyed, but those in LDCs and ISF would 
be unaffected. 

The over-production in this range need not be extreme 
as they survive in around 60 per cent of star forming regions 
(LDCs and ISF), but up to 50 per cent of those that form 
may be destroyed in HDCs. 

Thus, all modes of star formation combined must pro- 
duce about 10 per cent of G-dwarfs and 5 per cent of M- 
dwarfs with a companion between 50 and 1000 AU. 

• A similar fraction of 10 - 15 per cent of binaries are 
'soft-intermediate'. Those produced in HDCs will be almost 
all destroyed, many will in LDCs, but those in ISF will re- 
main. 

Following the above arguments, 40 per cent of star for- 
mation (HDCs) cannot produce soft-intermediates, in 40 per 
cent of star formation half of those that form may be de- 
stroyed (in LDCs), and in 20 per cent (ISF) all survive. 

Thus, all modes of star formation combined must pro- 
duce at least 10 per cent of G-dwarfs and 5 per cent of M- 
dwarfs with a companion between 10^ and 10* AU. Although 
the fraction formed could be significantly larger. 

• That leaves around 20 per cent of binaries that are 
'always soft' and cannot survive (or even be formed) in any 
cluster. At first inspection, it appears that such binaries 
must all form in ISFQ- 

If these binaries are produced by the ISF mode, 20 per 
cent of star formation must produce 20 per cent of the to- 
tal number of binaries. This implies that almost all isolated 
star formation must produce a binary with a separation 
> 10* AU. 

However, it is difficult to see how even isolated star 
formation can produce binaries with separations > 10* AU. 
Isolated star forming cores only have radii of ~ 0.1 pc out 
t o the point at which they m erge with the background (e.g. 
Ward- Thompson et al.|[2007l . and references therein) and so 
even a companion forming at the very limit of the core from 
the primary would only have a 2 x 10* AU separation (and 
surely disc fragmentation could not work at such distances) . 
Thus the origin of binaries > a few x 10"^ AU is a mystery. 

Clearly, the initial binary separation distribution can- 
not be identical to the field. However, the form and uni- 
versality (or otherwise) of the initial binary separation dis- 
tribution remains unclear. Note that there are two possible 
'initial' binary separation distributions (BSDs). Firstly, the 
'primordial' BSD produced as the immediate outcome of star 
formation, i.e. the distribution that emerges from Class O/I 
sources. Secondly, the 'initial' BSD that evolves rapidly from 
the primordial BSD due to circularisation and interactions 
with discs. It is the e volution from primordial to initial BSDs 
that eigenev o lution j[Kroupalll995b) attempts to capture. 

iKroupal (|l995al fa) described a potential initial BSD in 
which the closer binaries (< 50 AU) have a field- like distri- 
bution, and wider binaries over-produced by a factor of ~ 2 
out to 10* AU. However, as we have seen, the formation of 
binaries as wide as 10* AU is problematic at best in HDCs 
and all will be destroyed in HDCs. We will return to the 



° We are investigating the possibility that they may be formed 
during the violent dissolu tion of clusters after gas expulsion 
llKouwenhoven et al.ll2009bl . in prep). However, it is very unclear 
at the moment if this is possible. 
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problem of the form of the initial BSD and whether it is 
universal in a future paper. 



6 CONCLUSIONS 

We use Ai'-body simulations to dynamically evolve star clus- 
ters rich in binary systems to examine the effect of dynam- 
ical interactions on the initial binary population. Our main 
conclusions are: 

• Binary processing occurs within only a few initial cross- 
ing times, therefore it is the initial density of a cluster 
that is of importance in determining the binary population. 
Clusters older than a few (initial) crossing times will have 
reached a dynamical equilibrium with their binary popula- 
tion, dynamically young clusters will have not. 

• Binaries in clusters cannot form with the field separa- 
tion distribution. Binaries with separations > 10'' AU are 
too wide to form in any cluster, and binaries with separa- 
tions > 10^ AU are all destroyed in probably half of clusters. 

• Binaries in clusters cannot form with the field binary 
fraction for systems with separations in the range 50 to a 
few thousand AU. Many intermediate binaries are destroyed, 
therefore they must be overproduced relative to the levels 
observed in the field. 

• Binaries in clusters must form with roughly the cor- 
rect binary fraction and separation distribution for systems 
closer than '-^ 50 AU. Hard binaries are relatively unaffected 
by dynamical evolution; small amounts of destruction and 
modification do occur, but not enough to seriously affect 
this population. 

• The central region of the Orion Nebula Cluster was ini- 
tially ~ 6 times smaller (~ 200 times denser) than it is to- 
day. The binary populations were established and well-mixed 
at this time when the crossing time was significantly shorter. 
Therefore we expect there to be only small differences be- 
tween the binary popula tions in the inner an d outer regions 
of Orion as is observed (|Kohler et al.|[200^ . However, fur- 
ther observations are required to help distinguish between 
different models for the current and past state of the ONC. 

A comparison of the results of star formation to the 
field distribution at separations < 50 AU probably should 
be made, as this population should not normally be af- 
fected by dynamical evolution. However, it is not clear if 
low- and high-density environments and isolated star for- 
mation should always produce the same distribution(s). 

It appears that star formation should - for all masses, 
including M-dwarfs - produce more binaries with separa- 
tions in the range 10^ - 10* AU than are observed in the 
field as such binaries are highly susceptible to destruction 
in clusters (depending on density). But again, it is not clear 
if low- and high-density environments and isolated star for- 
mation should always produce the same distribution(s). 

However, star formation in clusters need not produce 
any binaries with separations > 10* AU, as such binaries - 
although observed in the field - just cannot survive in clus- 
ters, even if they can form. Their origin is a major problem 
in star formation and star clusters as it cannot be explained 
if most stars form in clusters. 
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